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I. CRITICAL COUPLING 

The critical coupling for the homogeneous case ilci is de- 
termined for the case of A^i = A^2 = N/2 with the linear 
density p = N/2L. The result (reproduced from the main 
text) is 



Mlcr = 



2(.9i2 - giig22)p 
ffii + 522 + 2512 



(1) 



For the harmonically trapped system, we use the Thomas- 
Fermi approximation for the ground state density conden- 
sate in the miscible regime when ri(0) > ilcr 



p{x) 



512+5 



(2) 



where we have used g = gn = 922 7^ 5i2- The chemical po- 
tential ^ is found by solving J p{x)dx = N. The spatial de- 
pendence of the density leads to the critical coupling strength 
ilcr quoted in the main text being spatially dependent 
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and the miscibility condition is fi(0) > fici (a^) for all x. We 
have numerically checked the validity of this expression for 
our trapping parameters within the local denxity approxima- 
tion. 



II. CONTROL PARAMETER 

The control parameter e [see Eq. (2) of the main text] for 
a harmonically trapped BEC is spatially dependent due to the 



inhomogeneity of the critical coupling in Eq. (O 
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where p{x) is given by Eq. (|2|. The phase transition now 
has a travelling front, beginning at a; = where the density 
is largest, and propagating outwards towards the edge of the 
BEC. The position of the front as a function of time is found 
from the solution of etr{x, t) = (2j] which yields 



XF{t) = ±y 
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It follows that the front moves with velocity 
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which diverges for \x-p\ — > 0. For very fast quenches, tq 0, 
the front is always faster than the speed of sound. For suffi- 
ciently long quenches there is a region X where the speed of 
the front is smaller than the local speed of sound 



Vs{x) 
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and domain formation is suppressed. The region X grows in 
size as rq increases, hence reducing the length of the conden- 
sate where domains form according to the homogeneous KZ 
mechanism. The increase of the size of the suppression region 
with the quenching time contributes to the modified scaling 
exponent for spatially inhomogeneous phase transitions. For 
a homogenous system with v-p <C Ws the scaling exponent for 
the binary BEC is predicted to be n = 4/3 lUt]. However, our 
result of ri k, 1/2 suggests that causaUty does not play a major 
role in the trapped system. 
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The imscibility-immiscibility phase transition in binary Bose-Einstein condensates (BECs) can be controlled 
by a coupling between the two components. Here we propose a new scheme that uses coupling-induced pattern 
formation to test the Kibble-Zurek mechanism (KZM) of topological-defect formation in a quantum phase 
transition. For a binary BEC in a ring trap we find that the number of domains forming the pattern scales as 
a function of the coupling quench rate with an exponent as predicted by the KZM. For a binary BEC in an 
elongated harmonic trap we find a different scaling law due to the transition being spatially inhomogeneous. We 
perform a "quantum simulation" of the harmonically trapped system in a ring trap to verify the scaling exponent. 



The formation of topological defects in symmetry-breaking 
phase transitions is a universal phenomenon relevant to many 
fields of physics, from cosmology [ij] to condensed matter 
0]. The Kibble-Zurek mechanism (KZM) jsUl^ is a theory 
relating the density of defects in the broken symmetry phase 
to the timescale of the transition. The scaling law predicted 
by the KZM has been demonstrated in simulations of phase 
transition dynamics ifisl - fisll . and has motivated many exper- 
iments HIt?, 19 -23il . However, the KZM scaling of defect 
density with the quench rate has not been verified in the labo- 
ratory. 

An experiment aiming to observe the Kibble-Zurek (KZ) 
scaling of defect production must have good control of the 
progress of the system through the phase transition. The ex- 
traordinary degree of flexibility and control available in ultra- 
cold gas experiments makes them leading candidates for the 
accurate testing of the KZM predictions. For example, the 
spontaneous formation of vortices in a Bose-Einstein con- 
densate (BEC) following the evaporative cooling of a single- 
component thermal Bose gas was suggested as a candidate 
system for studying the KZM ll24ll . Recently two experi- 
ments have observed this phenomenon, although they were 
not suited for testing the KZM j 25s, 26]. As Bose-Einstein 
condensation is a classical (thermal) phase transition, the for- 
mulation of a KZ scenario for this system requires the con- 
trol of thermodynamic parameters such as the temperature T 
and/or the chemical potential i^l. This is difficult for an iso- 
lated quantum system such as an ultra-cold gas, as the system 
thermalises under its own dynamics rather than through a cou- 
pling to a thermal reservoir 

The KZM, however, has been adapted to describe defect 
formation in quantum phase transitions ll27l - l29ll . in which a 
Hamiltonian parameter is ramped through a quantum critical 
point. Hamiltonian quenches are relatively straightforward to 
control in an isolated quantum system, and thus suggests that 
quantum phase transitions are strong candidates for quantita- 
tive KZM studies with ultra-cold gases. An important exper- 
iment on this topic was a study of the dynamics of a spin- 
1 BEC following a magnetic field quench js^l . resulting in 
the formation of topological defects such as spin-vortices and 
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polar-core vortices ll30ll . This experiment was analysed in the 
context of the KZM in lfl6l - [l8ll . but the KZ scenario in this 
system is complicated by the unclear role of dipolar interac- 
tions ifsill and the difficulties involved in imaging and count- 
ing the topological defects Isoll . A quantum quench from the 
Mott insulator to superfluid state in a optical lattice was re- 
cently reported in II32I1 . but did not demonstrate KZ scaling. 

In this Letter we formulate a straightforward, experimen- 
tally realistic Kibble-Zurek scenario for a coupling-induced 
miscibility-immiscibility quantum phase transition in a binary 
Bose-Einstein condensate jsll . The observable is the num- 
ber of domains formed in the immiscible phase. The domain- 
walls are stable and long-lived in a ring trap, and may be easily 
detected using absorption imaging Il34ll35ll . The transition is 
achieved by reducing the coupling between the two compo- 
nents provided by a microwave or laser field, ensuring precise 
control of the quench. 

We find that a coupling quench of a binary BEC in a ring 
trap confirms the KZ theory predictions for the resulting num- 
ber of defects. However, we find that the same transition in a 
quasi- ID harmonically trapped BEC yields a different scaling 
exponent for the defect density as compared to the ring BEC 
and the KZ prediction. We attribute the difference to the mo- 
tion and decay of the defects. To confirm this we design a spa- 
tially dependent quench of the Hamiltonian parameters in the 
ring trap such that the experiment has the same characteristics 
as the spatially inhomogeneous phase transition occurring in 
a harmonically trapped binary BEC. 

The physics of binary BECs has been extensively studied 
in previous work illSl- Here we consider a binary BEC 
formed by a single atomic species with two hyperfine states 
{i = 1,2) that are linearly coupled with time-dependent am- 
plitude ri(t), such as may be achieved with resonant two- 
photon microwave coupling ll4lll . We restrict ourselves to 
an elongated system where the spin dynamics are confined 
to the X dimension. We integrate out the transverse degrees of 
freedom (assumed to be harmonically trapped with frequency 
Lo±), and utilise the rotating-wave approximation (RWA) to 
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give the one-dimensional Hamiltonian 
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where Hq = — + Vix) is the single-particle Hamil- 

tonian, and tpi = ipi (x) is the Bose field operator for com- 
ponent i. The detuning of the transition is h6, which we 
subsequently set to zero. The nature of the ground state of 
the uncoupled system with n{t) = is determined by the 
parameter ll39ll A = .911322/312' where 312 (gu) is the ID 
inter-species (intra-species) interaction constant, with gij = 
2h'^aij /{maj_), a± — {h/mu}±y/^ and are the scatter- 
ing lengths. The Gross-Pitaevskii ground state ll42ll of Eq. ^ 
with V{x) — n{t) = is found to be miscible for A > 1 

Here we consider a system with A < 1 that is phase- 
separated for rt{t) = 0, but becomes miscible above a density 
dependent critical coupling strength ficr lUltl- If we begin an 
experiment in the ground state with f2(0) > ficr, the system 
will phase separate as il{t) is ramped to zero. If this occurs 
sufficiently quickly, the system cannot adiabatically follow its 
ground state, and a spatially random pattern of domains con- 
sisting entirely of atoms in either state 1 or 2 will form. Our 
proposal will allow for a test of the KZM by controlling the 
rate at which rt{t) is reduced to zero, and counting the final 
number of domains. 

We first consider a periodic uniform system = 0] of 

length L as realised by a binary BEC in a ring trap ll43l - l45ll . 
Within the mean-field approximation, the energy spectrum of 
the Hamiltonian ([T]i has agar 

in the miscible regime 11461 14711 . where hflcr ~ 2(9^3 ^ 
5ii522)p/(3ii + ff22 + 2^12) is the coupling strength that de- 
fines the quantum critical point and p is the linear atom density 
ll48ll . To formulate the Kibble-Zurek scenario for this system, 
we define a control parameter e{t) measuring the distance of 
the system from the critical point 



e{t) = 1 - n{t)/n^ 



(2) 



In the thermodynamic limit, when e{t) the correlation 
length ^ and the relaxation time t diverge as £^ ~ £^o/\e\'^ , t = 
To/jel"^, where the spatial and dynamical critical exponents 
(1/ and z respectively) depend only on the universality class 
of the transition. Due to the divergence of the relaxation time 
T as e 0, we expect the system evolution to become non- 
adiabatic when the time needed to adjust to an external change 
is equal to the timescale at which the coupling changes ifioll 



r(t) = e{i.)/e{i), 



(3) 



which defines the freezing time, i. Choosing a linear ramp 
il{t) = max [0, 2f2cr (1 — t/'^'Q)] for the coupling, where tq 
is the quench time, we find that the correlation length at the 
freezing time is ^(t) = Co('''q/''"o)''^'"'^^'^^''- For our system 
the mean-field critical exponents are v — 1/2 and z = 1 
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FIG. 1. (Color online) Examples of domain formation in the density 
of one hyperfine component for (a) a uniform system with tq = 
63 ms and TV = W atoms and (b) a trapped system with similar 
parameters. The red dashed vertical lines indicate the time where 
Q{t) — Q,ci, while the solid green vertical lines indicate the end of 
the quench 0,{t) — 0. The density of the second component (not 
shown) is complementary to the density of the first component. 



IHTIl . with ^0 = Cs/"\/2 where the spin healing length = 
h/y/2mpgs, and To = fi./2.gsp with 5^ = (.gii+522-2gi2)/2. 
The prediction for the mean number of domains at the end 
of the quench is 



= L/at) = (L/eo)(To/TQ)i/^ 



(4) 



To test this prediction, we perform quantum dynamical sim- 
ulations of the proposed experiment for a range of quench 
times TQ using the truncated Wigner method ll49ll . This is 
equivalent to adding a half quantum of noise per mode to the 
initial mean-field wave function, and then numerically solv- 
ing the time-dependent Gross-Pitaevskii equation for the bi- 
nary system IdOII . Quantum expectation values for the dy- 
namics then can be calculated by averaging over an ensem- 
ble of trajectories. We interpret each trajectory as an indi- 
vidual experimental realisation (4^, and count the number of 
domains at the end of the evolution. A typical trajectory is 
shown in Fig.lTJa), where domain formation is observed once 
il(t) < flcr- At times when < ^l{t) < flcr the domains 
can drift and merge ifsill . however once ^{t) =0 the pattern 
of defects is stable, allowing for their unambiguous counting 
B. 

For the ring system we choose parameters iV = 5 x 10''' 
atoms, L ~ 96 p,m, an = a22 = ai2/2 = 1.325 nm, and 
uj±/2Tr = 2 kHz. With such parameters we have = 1-16 
/im, ^0 = 0.82 /im and tq = 1.8 ms. The spin healing length 

is about two-thirds of the transverse system size, ensuring 
that the quasi- ID approximation is valid. We simulate quench 
times Tq over three decades in the range [0.1, 125] ms, and 
plot the mean number of domains versus the quench time 
Tq in Fig. |2] We fit a power law oc Tq" to the data for 
Tq > 2 ms, and find n ~ 0.341 ± 0.006 in good agreement 
with the KZM prediction of n = 1 /3. Thus a binary BEC in a 
ring trap is a strong candidate system for an experimental test 
of the KZM. 

We note that our numerical results deviate from the KZM 
prediction for rapid quenches with tq < 2 ms. For these 
quench times we find the number of domains is still decreas- 
ing at the end of our integration time (see inset of Fig. [2]) and 
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FIG. 2. (Color online) Mean number of domains formed A^d as a 
function of the quench time rq for a BEC in a ring trap. Blue circles 
are for a spatially homogeneous quench, and the red diamonds are for 
a spatially inhomogeneous quench simulating a harmonically trapped 
system. The open symbols average over 100 trajectories, while the 
closed symbols average over 1000 trajectories and are used for fitting 
the scaling exponent n. The error bars indicate the standard error. 
The linear fits to the data yield a scaling exponent n = 0.341 ± 
0.006 for the homogeneous quench, and n = 0.497 ± 0.015 for 
the inhomogeneous quench. Data for the inhomogeneous quench 
are shifted upwards by a multiplicative factor for clarity. The top 
right inset shows the same data on a linear scale. The bottom left 
inset shows the time evolution of the average number of domains for 
quenches with tq < 2 ms, for which the mean number of domains 
is still decreasing att — 600 ms. 



thus the mean number of domains plotted in Fig. |2] overes- 
timates the true final number. Numerical instabilities, com- 
bined with the known limitations of the truncated Wigner 
method, prevent us from further extending the integration time 
ll49ll . The upper limit to the number of domains is apparent, 
being iVI^^'' w « 80. 

We now perform simulations of a quasi- ID harmonically 
trapped system [^(x) = muj'^x'^ /2] with a;/27r = 5 Hz and 
other parameters as for the uniform ring trap. For a system 
with nonuniform density p{x) the critical coupling strength 
ilcr{x), and hence the control parameter etr{x,t), are spa- 
tially dependent and therefore the quantum phase transition 
is spatially inhomogeneous. The denser central region of the 
condensate enters the immiscible phase earlier, and a mo ving 
front separates it from the parts still in the miscible phase 15211 
as seen in Fig.fljb). 

The stochastic nature of the simulations means that it is dif- 
ficult to distinguish domains from density variations near the 
edge of the condensate, and hence the counting of domains 
after the quench in the trapped system is problematic. To 
overcome this issue, we restrict our counting to the central re- 
gion, where the initial density is greater than a threshold value 
Pcut = 7P(0). Figure [3] shows results for the final number of 
domains as a function of quench time for a range of 7. The 
inset of Fig. [3] shows the scaling exponent resulting from the 
fitting of these curves to A^d oc Tq " as a function of 7. The 
approximately constant value of n w 0.47 for 7 > 0.3 sug- 
gests that variation of the scaling exponent with 7 is mostly 



10 



o Y=0.1 [R = 0.g5R.^p] 
y= 0.2 [R = 0.90 R.^p] 
^ 7 = 0.4[R = 0.77R.^p] 
T=0.6 [R = 0.63R.^p] 




10*^ 



10 



Xq (ms) 



10^ 



FIG. 3. (Color online) Mean number of domains formed A^d ver- 
sus quench time rq for a harmonically trapped BEC. The different 
colors (symbols) correspond to different domain-counting threshold 
densities 7 — pcut/p(0) (see text). Filled symbols represent ensem- 
bles of 1000 trajectories used for estimating the scaling exponent n. 
The left inset shows the value of the exponent n as a function of the 
density threshold 7. The dashed (dot-dashed) line indicates the value 
of n obtained from a simulation of a homogeneous (inhomogeneous) 
coupling quench in a ring trap as in Fig. [2] (see text). The shaded 
region in the right inset illustrates the counting region in the BEC 
density. The flattening of the domain number at long quench times 
for 7 = 0.1 and 7 — 0.2 is due to the unavoidable miscounting of 
density fluctuations as domains. 



due to the miscounting of domains. 

To help understand the inhomogeneous phase transition, 
we model a "quantum simulation" of the spatially uni- 
form quench in a harmonically trapped BEC by implement- 
ing a spatially dependent ramp of the coupling strength 
in a ring BEC. This is achieved by choosing U{x,t) = 
max[0, 2f]cr(l - t/TQ)p{Q)/ p{x)] with g = gn = 322 
312/2, where p{x) is the ground state density of the harmon- 
ically trapped system we are modelling. For such a quench 
the control parameter in the ring is spatially dependent and 
is identical to eti{x, t) for the harmonically trapped BEC. To 
avoid the divergence of €tr{x, t) where p{x) — s- 0, we choose 
the circumference of the ring such that L < 2Rtf, where 
i?TF is the Thomas-Fermi radius of the trapped system being 
simulated ifsill . We plot the number of domains versus quench 
time in Fig. |2l and find clear evidence for power-law scaling 
with an exponent n = 0.497 ± 0.015. This "quantum simula- 
tion" demonstrates that our result for the scaling exponent in 
the harmonic trap is robust, but we are unable to obtain this 
result analytically ifsill . 

We identify two effects potentially contributing to the in- 
crease in the scaling exponent in the harmonic trap compared 
to the ring trap. First, during an inhomogeneous phase transi- 
tion the difference between the coupling and the critical cou- 
pling strength is spatially dependent, and this introduces a 
preferred direction for the movement of domains. The break- 
ing of translation invariance leads to a larger annihilation rate 
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of domains than in a homogeneous phase transition. During 
longer quenches domains have more time to annihilate or es- 
cape the counting region, resulting in an increase of the ob- 
served scaling exponent 115 ill . 

Second, it has been proposed that when the front velocity 
vf{x) is less than the local speed of sound Vs{x), the parts 
of the system in the symmetry -broken phase can influence the 
symmetry breaking in the region undergoing the transition, 
leading to a suppression of domain formation. This was previ- 
ously discussed with regard to soliton formation in a thermally 
quenched Bose gas in Ref. iIstIi . Below a certain tq there is 
a spatial region for which the velocity of the front is less than 
the speed of sound, and this increases in size with increasing 
TQ. Thus slower quenches result in a larger spatial region in 
which domain formation is suppressed, and this contributes to 
a larger scaling exponent. 

Finally, we consider the feasibility of the experiments we 
propose. No pair of hyperfine states of ^^Rb or ^■^Na natu- 
rally satisfies the immiscibility criterion with a spin healing 
length sufficiently small to allow for the formation of mul- 
tiple domains. However, the two hyperfine states |1) = = 
l,mF = +1) and |2) ^ \F ^ 2,mF = -1) of ^'^Rb have 
.911 ~ .922, and exhibit an interspecies Feshbach resonance 
that could be utilised to tune gi2 llssi Isill . We estimate that 
it is possible to attain A w 0.8 while keeping inelastic losses 
sufficiently small to allow enough time to realise the proposed 
experiment. For a ring BEC of 5000 atoms with uj ±/2Tr = 2 
kHz and a circumference of 50 /im, A = 0.8 implies = 0.4 
Aim, To = 351 ms and iV'"^'^ w 50. 

The miscible-immiscible phase transition in binary conden- 



sates can also be controlled via the spin-orbit coupling of neu- 
tral atoms 1I55I1 . In this situation the phase transition to the 
immiscible state is achieved by ramping up the intensity of 
two slightly detuned lasers coupling two hyperfine levels of 
= 1) of **^Rb. This scheme has the significant advantage 
of being able to reach the strongly immiscible phase without 
suffering from the inelastic atom losses common near a Fesh- 
bach resonance ll56ll . 

For both schemes the stability of domains is ensured far 
from the transition, as each component acts as an effective 
potential for the other. In the strongly immiscible regime 
A < 0.95 the kinetic and thermal energy are not sufficient 
to overcome the barriers provided by the domain pattern. 

In conclusion, we have shown that the number of domains 
arising in coupling-induced pattern formation in a binary BEC 
in a ring trap scales as predicted by the Kibble-Zurek mech- 
anism. Recent demonstrations of ring BECs ll43l - l45ll . com- 
bined with the experimental feasibility of the scheme, make 
it an excellent candidate for testing the Kibble-Zurek theory. 
We have also verified that a scaling law exists for harmoni- 
cally trapped BECs, allowing for a qualitative test of the KZM 
in this system. 
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